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' A compelling framework to explain dark matter and electroweak baryogenesis is supersymmetry with light scalar 

^Nj , top quarks (stops) and a small mass difference between the stop and the lightest neutralino. In this work, the stop 

detection capability at the ILC for small mass differences between the stop and the neutralino is studied. The analysis 
r-{ I is based on a fast and realistic detector simulation. Significant sensitivity for mass differences down to 5 GeV is 

obtained. It is discussed how the relevant parameters of the scalar tops can be extracted and used to compute the 
dark matter density in the universe. 
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1. INTRODUCTION 

> 

' The origin and stabihzation of electroweak symmetry, the nature of dark matter and the generation of the baryon 
QQ ' asymmetry (baryogenesis) in the universe suggest the existence of new symmetries within the reach of next generation 
colliders. A particular attractive extension of the Standard Model that addresses all these issues is supersymmetry 
with TeV-scale supersymmetry breaking. 

A long history of experimental observations has corroborated the evidence for dark matter in the universe, cul- 
O ^ minating in the recent accurate determination by the WMAP satellite, in combination with the Sloan Digital Sky 
' ; Survey (SDSS) Q, ^cdm/i' = OAUGt^ofiH at the 95% C.L. Here ricDM is the dark matter energy density normal- 
' ized to the critical density and h is the Hubble parameter in units of 100 km/s/Mpc. Supersymmetry with i?-parity 
conservation provides a natural dark matter candidate, which in most scenarios is the lightest neutralino. 
^ ' Electroweak baryogenesis is based on the concept that the baryon asymmetry is generated at the electroweak 
. phase transition. While in the Standard Model the phase transition is not sufficiently strongly first order and there 
^ ' is not enough CP violation, supersymmetry can alleviate both shortcomings. A strong first-order phase transition 
- - ' can be induced by loop effects of light scalar top quarks (stops) to the finite temperature Higgs potential. In much of 
the parameter space of interest for electroweak baryogenesis, the light stop is only slightly heavier than the lightest 
neutralino, thus implying that stop-neutralino co-annihilation is significant. The co-annihilation contribution is 
important to meet the WMAP experimental results shown above. In the co-annihilation region, the stop-neutralino 
mass difference is smaller than 30 GeV Q , making a discovery of the stops at hadron colliders difficult Q . 

This work investigates the capabilities of a future international e+e" hnear collider (ILC) to discover light stops 
and measure their properties. This study extends the reach of an earlier work 4] to the region of small stop-neutralino 
mass differences. The experimental analysis for the ILC is used to assess possible conclusions about the nature of 
dark matter. 



2. EXPERIMENTAL SIMULATION 

In this work, the production of light stops at a 500 GeV linear collider is analyzed, using high luminosity C ~ 
500 fb~^ and polarization of both beams. For small mass differences Am — TOj^ — m-^o, a light stop decays dominantly 
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Table I: Cross-sections for the stop signal and Standard Model background processes for ^/s = 500 GeV and different polar- 
ization combinations. The signal is given for the stop mixing angle cosSj = 0.5. Negative/positive polarization values refer to 
left-/right-handed polarization, respectively. 



Process 


Cross-section [pb] 


P(e-)/P(e+) 


0/0 


-80%/+60% 


+80%/-60% 


iiil mi^ = 


120 GeV 


0.115 


0.153 


0.187 


'Tlti = 


140 GeV 


0.093 


0.124 


0.151 


"^ti = 


180 GeV 


0.049 


0.065 


0.079 




220 GeV 


0.015 


0.021 


0.026 


w+w- 




8.55 


24.54 


0.77 


zz 




0.49 


1.02 


0.44 


Wev 




6.14 


10.57 


1.82 


eeZ 




7.51 


8.49 


6.23 


qq, q^t 




13.14 


25.35 


14.85 


ti 




0.55 


1.13 


0.50 


2- phot on, pt > 5 GeV 


936 


936 


936 



into a charm quark and the lightest neutralino, ti — > cxi- The signature for stop pair production at an e+e collider, 

e+e- ^hil^ccx^.Xi, (1) 

is two charm jets plus missing energy. For small Am, the jets are relatively soft and separation from backgrounds 
is very challenging. Backgrounds arising from various Standard Model processes can have cross-sections that are 
several orders of magnitude larger than the signal, so that even small jet energy smearing effects can be important. 
The possibility of efficient charm tagging plays an important role in the signal identification and depends on the 
vertex detector performance. Thus it is necessary to study this process with a realistic detector simulation. Signal 
and background events are generated with Pythia 6.129 0] with CiRCE for realistic beamstrahlung together 
with additional private code for the stop signal generation previously used in Ref. P|. The detector simulation is 
based on the fast simulation Simdet 0, describing a typical ILC detector. 

Tab. m lists the cross-sections for the signal process and the relevant backgrounds. They have been computed 
with the Monte-Carlo code used in Ref. jj and by Grace 2.0 0, with cross-checks to CompHep 4.4 jl3| where 
applicable. A minimal transverse momentum cut, pt > 5 GeV, is applied for the two-photon background, to avoid 
the infrared divergence. 

In the first step of the event selection, the following cuts are applied: 

4 < A^chargod tracks < 50, Pt > 5 GcV, 

I C0s6'thrust| < 0.8, blong,tot/Ptot| < 0.9, (2) 

-Bvis < 0.75^5, minv < 200 GeV. 

The cut on the number of charged tracks removes most leptonic background and part of the ti background. By requir- 
ing a minimal transverse momentum pt, the two-photon background and back-to-back processes like qq are largely 
reduced. The signal is characterized by large missing energy and transverse momentum from the two neutralinos, 
whereas for most backgrounds the missing momentum occurs from particles lost in the beam pipe. Therefore, cuts 
on the thrust angle ^thrust, the longitudinal momentum piong.tot, the visible energy E^is and the total invariant mass 
TOinv are effective on all backgrounds. As can be seen in Tab. ^ the various background are substantially reduced 
after these preselection cuts, while roughly around 70% of the signal is preserved. 

After generating large event samples with the preselection cuts for the various backgrounds as listed in Tab. ^ 
the following final event selection cuts are applied to further improve the signal-to-background ratio: 

1. Number of jets Ajcts = 2. Jets are reconstructed with the Durham algorithm with the jet resolution parameter 
2/cut — 0.003 X ^/s/ E^is- The cut reduces substantially the number of W and quark pair events. 
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Table II: Background event numbers and titl signal efficiencies in % (for various mj^ and Am in GeV) after preselection and 
each of the final selection cuts. In the last column the expected event number are scaled to a luminosity of 500 fb~^. The cuts 
are explained in the text. 









After 














Scaled to 


Process 




Total 


presel. 


cut 1 


cut 2 


cut 3 cut 4 


cut 5 


cut 6 


500 fb"^ 






210,000 


2814 


827 


28 


25 


14 


14 


8 


145 




zz 


30,000 


2681 


1987 


170 


154 


108 


108 


35 


257 




Weu 


210,000 


53314 


38616 


4548 


3787 


1763 


1743 


345 


5044 




eeZ 


210,000 


51 


24 


20 


11 


6 


3 


2 


36 




qq,q^t 


350,000 


341 


51 


32 


19 


13 


10 


8 


160 




ti 


180,000 


2163 


72 


40 


32 


26 


26 


25 


38 




2-photon 


8 X 10*^ 


4061 


3125 


3096 


OOO 


A no 








< 164 


m;^ = 140 


Am = 20 


50,000 


68.5 


48.8 


42.1 


33.4 


27.9 


27.3 


20.9 


9720 




Am = 40 


50,000 


71.8 


47.0 


40.2 


30.3 


24.5 


24.4 


10.1 


4700 




Am = 80 


50,000 


51.8 


34.0 


23.6 


20.1 


16.4 


16.4 


10.4 


4840 


mi^ = 180 


Am = 20 


25,000 


68.0 


51.4 


49.4 


42.4 


36.5 


34.9 


28.4 


6960 




Am = 40 


25,000 


72.7 


50.7 


42.4 


35.5 


28.5 


28.4 


20.1 


4925 




Am = 80 


25,000 


63.3 


43.0 


33.4 


29.6 


23.9 


23.9 


15.0 


3675 


mi^ = 220 


Am = 20 


10,000 


66.2 


53.5 


53.5 


48.5 


42.8 


39.9 


34.6 


2600 




Am = 40 


10,000 


72.5 


55.3 


47.0 


42.9 


34.3 


34.2 


24.2 


1815 




Am = 80 


10,000 


73.1 


51.6 


42.7 


37.9 


30.3 


30.3 


18.8 


1410 



2. Large missing energy, E^is < 0■4:^/s. This is effective against W^W^ , ZZ and di-quark events. In addition, 
a window for the invariant jet mass around the W-boson mass, 70 GcV < mjot.inv < 90 GeV, is excluded to 
reduce the large Weiy background. 

3. qq events are removed by requiring a minimal acoUinearity angle cos(/)aco > —0.9. 

4. Cutting on the thrust angle, | cos6'thiust < 0.7|, reduces W boson background. 

5. A strong cut on the transverse momentum, pt > 12 GeV, completely removes the remaining two-photon events. 

6. The largest remaining background is from e+e" — > Wev. It resembles the signal closely in most distributions, 
e.g. as a function of the visible energy, thrust or acoUinearity. Only by increasing the invariant jet mass window 
from cut 2 to (60 GeV < mjet.mv < 90 GeV), the signal-to-background ratio is improved, but at the cost of 
losing a substantial amount of the signal. In addition, the signal selection can be enhanced by charm tagging, 
which is implemented based on the neural network analysis described in Ref. [3 ■ The neural network has been 
optimized to reduce the Weu background while preserving the stop signal for small mass differences. 

The resulting event numbers, scaled to a luminosity of 500 fb"^, and the signal efficiencies are listed in Tab. ITU 
After the final selection, the tit\ signal event numbers are of same order as the remaining background, N ^ 0(10'*). 

To explore the reach for very small mass differences Am = — to-^o , signal event samples have been generated 
also for Am = 10 GeV and 5 GeV, see Tab. IIIII The signal efficiency drastically drops for Am — 5 GeV, as a result 
of the Pt cut (cut 5). An optimization of the event selection for very small Am will be addressed in future work. 

Based on the above results from the experimental simulations, the discovery reach of a 500 GeV e+e~ collider can 
be estimated, see Fig. ^ To cover the whole parameter region, the signal efficiencies for the parameters points in 
Tab, mil are interpolated. Then, the signal rates S are computed by multiplying the efficiency e obtained from the 
simulations with the production cross-section for each point (m^^ , m^o ) . Together with the number of background 
events B, this yields the significance S/ VS + B. The green area in the figure corresponds to the 5a discovery 
region, S/\/S + B > 5. As evident from the figure, the ILC can find light stop quarks for mass differences down to 
Am ~ 0(5 GeV), covering the complete stop-neutralino co-annihilation region and beyond. 
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Table III: Signal efficiencies for titl production after final event selection for different combinations of the stop mass m^^ and 

mass difference Am — m,- — m-o . 

1 



Am 


mi^ = 120 GeV 


140 GeV 


180 GeV 220 GeV 


80 GeV 




10% 


15% 


19% 


40 GeV 




10% 


20% 


24% 


20 GeV 


17% 


21% 


28% 


35% 


10 GeV 


19% 


20% 


19% 


35% 


5 GeV 


2.5% 


1.1% 


0.3% 


0.1% 




Figure 1: Discovery reach of linear collider with 500 fb~^ luminosity at — 500 GeV for production of light stop quarks, 
e^e~ ^ iiii ^ ccxiXi- The results are given in the stop vs. neutralino mass plane. In the dark shaded region, a 5a 
discovery is possible. The region where m^o > mj^ is inconsistent with a neutralino LSP, while for mj^ > Mw + rrih + Xi the 
three-body decay ii W^bxi becomes accessible and dominant. In the light shaded corner to the lower left, the decay of 
the top quark into a light stop and neutralino is open. Also shown are the parameter region excluded by LEP searches [li| 
(white area in the lower left) and the Tevatron light stop reach 0| (dotted lines) for various integrated luminosities. 



3. STOP PARAMETER DETERMINATION 

The discovery of light stops would hint towards the possibility of electroweak baryogenesis and may allow the co- 
annihilation mechanism to be effective. In order to confirm this idea, the relevant supersymmetry parameters need 
to be measured accurately. In this section, the experimental determination of the stop parameters will be discussed. 

For definiteness, a specific MSSM parameter point is chosen: 

= -992 GeV^ At = -1050 GeV, Ah = 112.6 GeV, = 320 GeV, 
m^^ = 4200 GeV, tan^ =5, M2 = 225 GeV, (j),, =0.2. 

The chosen parameters are compatible with the mechanism of electroweak baryogenesis, generating the baryon 
asymmetry through the phase of /i. They correspond to a value for the dark matter relic abundance within the 
WMAP bounds, ^cuMh^ = 0.1122. The relic dark matter density has been computed with the code used in 
Ref. In this scenario, the stop and lightest neutralino masses are m^^ = 122.5 GeV and m^o = 107.2 GeV, and 
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the stop mixing angle is cos6'i — 0.0105, i.e. the hght stop is almost completely right-chiral. The mass difference 
Am = mi^ — m^o = 15.2 GeV lies within the sensitivity range of the ILC. 

The measurement of iiil production cross-section for different beam polarization allows to extract both the mass 
of the hght stop and the stop mixing angle 0|. Here is it assumed that 250 fb~^ is spent each for P(e~)/P(e+) = 
— 80%/+60% and +80%/— 60%, where negative/positive polarization degrees indicate left-/right-handed polarization. 
Besides the statistical errors, the following systematic errors are taken into account: mass measurement of the lightest 
neutralino ^rrij^o = 0.1 GeV (see Ref. for details); determination of the polarization degree (5P(e^)/P(e^) = 0.5%; 
error in the integrated luminosity 5C/C = 5 x 10""*; theoretical uncertainty for background simulation 5B/B = 
0.3%; stop hadronization and fragmentation 0.5-1.0%; charm fragmentation and tagging efficiency ^ 0.5%; detector 
calibration ~ 0.5%; beamstrahlung uncertainty from [l6j. 

Each of the two cross-section measurements for P(e~)/P(e+) = — 80%/-|-60% and +80%/— 60% correspond to 
a band in the parameter plane of the stop mass and mixing angle, see Fig. El Combining the two cross-section 
measurements, the stop parameter are determined to 

mi^ = (122.5 ± 1.0) GeV, cos 6*4 < 0.074 ^ sin6't > 0.9972. (4) 



4. CONCLUSIONS 

Using the data from the experimental stop analysis together with estimated errors for measurements in the neu- 
tralino/chargino sector, see Ref. |l5j |. the expected cosmological dark matter relic density can be computed. The 
mass of the heavier stop <2 is too large to be measured directly, but it is assumed that a limit of m^^ > 1000 GeV 
can be set from collider searches. 

All experimental errors are propagated and correlations are taken into account by means of a analysis. The 
result of a scan over 100000 random points for the scenario eq. ((SJ is shown in Fig. 13 The horizontal bands depict the 
relic density as measured by WMAP |1||, which is at la level 0.104 < ricDM^^ < 0.121. The collider measurements 
of the stop and chargino/neutralino parameters constrain the relic density to 0.086 < i^cDM^^ < 0.143 at the Itr 
level, with an overall precision of the same order of magnitude, though somewhat worse, than the direct WMAP 
determination. The uncertainty in the computation of the dark matter relic density from ILC precision measurements 
is dominated by the measuremenst of the lightest stop mass. 
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Figure 3: Computation of dark matter relic abundance JIcdm/i^ taking into account estimated experimental errors for stop, 
chargino, neutralino sector measurements at the ILC. The black dots correspond to a scan over the la (Ax^ < 1) region 
allowed by the experimental errors, as a function of the measured stop mass, with the red star indicating the best- fit point. 
The horizontal shaded bands show the la and 2a constraints on the relic density measured by WMAP. 
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